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Abstract. A numerical approach is employed to explain transport characteristics 
in realistic, quantum Hall based Aharonov-Bohm interferometers. First, the spatial 
distribution of incompressible strips, and thus the current channels, are obtained 
applying a self-consistent Thomas-Fermi method to a realistic heterostructure under 
quantized Hall conditions. Second, the time-dependent Schrodinger equation is solved 
for electrons injected in the current channels. Distinctive Aharonov-Bohm oscillations 
are found as a function of the magnetic flux. The oscillation amplitude strongly 
depends on the mutual distance between the transport channels and on their width. At 
an optimal distance the amplitude and thus the interchannel transport is maximized, 
which determines the maximum visibility condition. On the other hand, the transport 
is fully suppressed at magnetic fields corresponding to half-integer flux quanta. 
The results confirm the applicability of realistic Aharonov-Bohm interferometers as 
controllable current switches. 
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1. Introduction 

The Aharonov-Bohm (AB) effect pQ is among the most significant and useful phenomena 
in quantum mechanics. The AB effect manifests itself in the interaction between a 
charged particle and an electromagnetic field, even if the local magnetic and electric fields 
are zero in that region. The necessary information is included in the vector potential 
A, which induces a phase shift in the wave function of the electron traveling along a 
specific path. In a double-slit system, or in a quantum ring (see Ref. [2] and references 
therein), the relative phase shift between two electrons traveling along different paths is 
A0 = 27r$/$o, where $ is the total magnetic flux enclosed by the path and $o — h/e 
is the magnetic flux quantum. The resulting current (and conductance) of the quantum 
ring is then a periodic function of $/$o- 

Recent low-temperature transport experiments [31 HI El El IE] performed at two- 
dimensional (2D) electron systems (2DESs) utilize the quantum Hall (QH) effect to 
investigate and control of the electron dynamics via their AB phase. An interesting 
difference between the original AB experiments and QH interferometers is the fact that 
in the latter the electron path itself may depend on the magnetic field B. To describe 
electron transport in QH interferometers, the single-particle edge-state approach [9] is 
common, but it neglects the dependency of the area enclosed by the current- carrying 
channels on the magnetic field [TO], as well as on the channel widths. However, as 
shown explicitly below, the actual paths can be obtained considering the full many-body 
electrostatics, which yields the spatial distribution of compressible and incompressible 
strips [TT] . 

The essential features in the observed AB oscillations in QH interferometers have 
been explained using edge-channel simulations and Coulomb interactions at the classical 
(Hartree) level [U [121 [13] • However, a complete theoretical picture of the observed 
phenomenon is still missing [HI [H]. To attain this, it would be particularly important 
to (i) describe the full electrostatics by handling the crystal growth parameters and the 
"edge" definition of the interferometer, and to (ii) supply this scheme with a dynamical 
study on electronic transport in the 2DES. 

The objective of this work is to take important steps towards comprehensive 
explanation of the AB characteristics in QH interferometers. First, we apply the three- 
dimensional (3D) Poisson equation and the Thomas-Fermi approximation for the given 
heterostructure [To] , taking into account the lithographically defined surface patterns. In 
this way we obtain the electron and potential distributions under QH conditions [HI [T7] . 
For completeness, we utilize this scheme for the real experimental geometry resulting 
from the trench-gating technique. Second, we determine a model potential describing the 
current channels, and use a time-dependent propagation scheme to monitor transport 
of a wave packet injected in the channel. We find distinct AB oscillations, whose 
characteristics dramatically depend on the channel widths and their mutual distance. In 
particular, we show that there is an optimal way to manipulate the visibility in realistic 
AB interferometers. 
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Figure 1. (Color online) (a) Schematic picture of the Hall resistance as a function of 
B. (b)-(d) Corresponding potentials (black solid lines), Landau levels (dashed lines), 
and wave functions (red solid lines). Here /i (gray solid lines) is the Fermi energy at 
equilibrium. The ellipses in (a) indicate the interval in B where the incompressible 
strips become well-developed (b), washed out (c), or leaky (d). The actual calculations 
considering realistic structures can be found in Refs. [23 [221 US] 



2. Device and electrostatics 

The current- carrying states in a QH device result from the Landau-level quantization 
followed by level bending at the edges where the Fermi energy crosses the levels. 
Thus, the transport takes place through the edge states. However, there has been 
substantial debate in the literature whether the current flows through the compressible 
or incompressible strips. Although the ballistic ID picture [18j [9], later attributed 
to compressible strips [19], applies well to the integer quantum Hall effect (IQHE), it 
requires bulk (localized) states [20] to explain the transitions between the QH plateaus. 
In contrast, the screening theory assumes that the current is carried by the scattering- 
free incompressible strips only if the widths of these strips (channels) are wider than 
the quantum mechanical length scales [21J. 

A schematic presentation of the Landau levels across a QH bar is given in Fig. [1] 
together with the Hall resistance. There is a steep potential variation in incompressible 
strips, where the Fermi energy falls in between two consecutive Landau levels. If the 
incompressible strip is well-developed and larger than the wave extent, the quantized 
Hall effect is observed [Fig. G^b)]. However, if the strip width Wis becomes comparable 
with the wave extent (or the magnetic length {t = ^Ti/eB) for the lowest Landau 
level) [Fig. [He)] the strip loses its incompressibility, and partitioning between channels 
becomes possible. Thus, it is possible to observe interference. This case is called as 
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the "leaky" incompressible state, which can be accurately determined by taking into 
account the collision broadened Landau level width and the extent of the wave function. 
Once the strip becomes even smaller than the wave extent, the classical Hall effect is 
observed [22] as shown in Fig. Q^d). Altogether, interference can be observed at the 
lower end of the quantized Hall plateau. 

We calculate of the electron density and the electrostatic potential at the layer of 
the 2DES self-consistently [23] by using the structural information from Goldman [I]. 
The dopant density, location of the interface for the 2DES, and the dielectric constant 
k (= 12.4 for GaAs/AlGaAs) are used as the input to calculate the total potential from 
V(r) = Konf (r)+Vi n t(r), where the confinement potential V con {(r = (x, y, z)) is composed 
of (i) the potential generated by surface pattern (corresponding to trench-gating [17]). 
(ii) donors, and (iii) surface charges. The interaction potential is calculated from the 
electron density n e \ by solving the Poisson equation, 



where the kernel K(f; f) takes into account the imposed boundary conditions. The 
explicit kernels considering different boundary conditions can be found in Refs. [2U OS] 
and references therein. The electron density is calculated within the Thomas-Fermi 
approximation from 



where D(E;r) is the density of states, and the Fermi occupation function 
f(V(r),E,T,fjL) depends on the chemical potential \i and temperature T. The above 
equations are solved self-consistently on a uniform 3D grid with open boundary 
conditions. 

It is important to note the advantages and constraints of the applied 
approximations. As a general note, the number of electrons is very large (of the order of 
10 3 ), so that semiclassical or density-functional-type approximations become reasonable 
in order to obtain a qualitative estimation for the positions and the widths of the current 
carrying channels. Their exact properties require formidable numerical approaches that 
are beyond the scope of this work. 

The magnetic field effects are included in our static calculations via the 
approximation of Landau-quantized density of states, although the wave functions 
should also depend on the position according to the total potential. However, this 
approximation can be justified for different regimes of compressibility. First, in 
compressible systems the external potential is well screened (for instance, similar to 
the potential at the center of Fig. QJd), so that the total potential is constant. Hence, 
the Landau wave functions are not modified at all, but the energy eigenvalues are 
shifted. Secondly, in incompressible systems the total potential varies linearly, and 
hence the wave functions preserve the Landau form; only their center coordinates are 
shifted slightly (less than the magnetic length). At the transition regions the potential 
may vary with higher-order corrections. However, the next (quadratic) correction also 
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preserves the Landau form, whereas the center coordinate remains unaffected. In total, 
the corrections to the Landau wavefunction are thus negligible up to the third order, 
which justifies the practical validity of our approximation in determining the electron 
density. To compensate for the suppression of the tunneling effects due to the Thomas- 
Fermi approximation we take into account the finite widths of the Landau wave functions 
as discussed below. 

The exact form of D(E) is essentially determined by the properties of the scattering 
mechanisms in the presence of strong disorder. These constraints are lifted by the 
experimental conditions, so that the number of electrons within the interferometer is 
more than 10 3 and the mobility of the samples is extremely high (typically, \i > 5 x 10 6 
cm 2 /Vs). Hence, in such samples the Landau level broadening can be neglected: 
T/(hu c ) = r / \eB I m*) <C 1, where T is the level width and hu c is the cyclotron energy. 
Otherwise, the level broadening has been proven to have an important influence on the 
(in) compressibility of the system [23] . 

We point out that the statistical description of a many-particle system by the 
grand-canonical ensemble assumes that the system is in contact with a reservoir, which 
can exchange particles with the system and determines the chemical potential. If the 
system is inhomogeneous, one usually considers a constant electrochemical potential in 
thermal equilibrium, which, however, is position dependent in the case of an imposed 
current, i.e., in a (global) non-equilibrium. 

To proceed with the calculation to the desired values of B and T, we impose 
periodic boundary conditions and replace the constant density of states in the absence 
of magnetic fields with a Gaussian-broadened one [16] that takes the quantization due 
to the magnetic field into account, together with disorder. In this procedure, we 
first increase T and smear the quantization effects with the Fermi function. Then, 
we set the desired value for B and decrease T stepwise to its target value. In each 
iteration step a relative accuracy threshold of < 10~ 6 is obtained for the density. 
This standard calculation scheme has been previously used to describe similar systems 
successfully (2D ESI QE1 E2] ■ 

In Fig. [2J we show the spatial distribution of incompressible strips at different 
magnetic fields. The sequence of figures at B = 7.6 . . . 8.8 T shows that the distance 
between the incompressible strips decreases, and their width increases, as a function of 
B. It can be expected that at low B when the strips are far apart, the AB oscillations 
are not present or they are weak, as the wave packet remains in the incoming channel. 
On the other hand, the strong overlap between the strips at B — 8.8 T suggests 
suppression as well - now due to the complete transport of the wave packet to the other 
channel. Interference would require clean splitting (partitioning) of the wave function 
between the channels. To confirm this phenomenological suggestion, and to study the 
effects of the channel distance and width on the conductance, we focus now on realistic 
modeling of current channels along the incompressible strips. This is followed by a 
dynamical study on the electron transport in the device. As the core outcome we can 
qualitatively describe the interference intervals depending on the system parameters 
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Figure 2. Spatial distribution of the incompressible strips (black) calculated at (a) 
B = 7.6 T, (b) 8.0 T, (c) 8.4 T, and (d) 8.8 T at T = 1 K. It is expected that only 
at B — 8.0 T one can observe Aharonov-Bohm oscillations, whereas in the other cases 
the scattering between the edge-states is prevented (see text). 

such as temperature, sample geometry, and heterostructure properties. 

Before introducing the model for dynamics, a brief discussion about the direction of 
the current along the incompressible strips is in order. It is well known that the current- 
carrying states in a QH device result from the Landau-level quantization followed by level 
bending at the edges where the Fermi energy crosses the levels. Thus, the equilibrium 
transport takes place through the chiral edge states, i.e., with different current 
directions at opposite edges, and the net current is zero. In transport experiments 
(including interference experiments), however, an external current is imposed and a 
Hall (electrochemical) potential difference - with the same slope - develops between 
two opposing edges [25]. This potential drop at the incompressible strips, that confines 
the external current in those regions, has been shown in several experiments [26ll27l i28j. 
In the following we utilize this picture in the transport calculations, so that the electrons 
flow to the same direction along the opposite channels. 

3. Electron transport 

The channels are modeled by a 2D potential profile consisting of two curved pipes 
following the shape of the incompressible strips in Fig. [2j The potential minima of the 
channels are shown as solid lines in Fig. |3j The distance between the left and right 
parts is varied such that at the two encountering points (bottom and up) the distance 
between the potential minima is d = . . . 1 (in atomic units). It is important to note 
that the channels are genuinely 2D according to the real device, and their cross sections 
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Figure 3. Snapshots of the electron density in a model for the Aharonov-Bohm 
interferometer. The two lowermost rows correspond to zero flux and half a flux 
quantum, respectively. 



2/2 

have a Gaussian form V croS s ~ ~ Vo &~ s > where s is the coordinate perpendicular 
to the channel, Vo = 20 is the channel depth, and c = 0.2 is the width parameter 
which is varied. The choice of a Gaussian profile is justified by considering the magnetic 
(parabolic) confinement, which is close to a Gaussian form at the bottom of the channel. 
On the other hand, in the upper part of the channel the selected profile allows "leaking" 
of the electron flow according to the experiment (see above). 

As an initial condition we set a single-electron wave packet in the lower-left corner 
of the device and, in order to emulate a source-drain voltage, we accelerate the wave 
packet using a linear potential with slope V\ = — 0.2Vo along the channel, which leads 
to an initial velocity of ~ 3.3 a.u. We monitor the density and the current density 
during the time propagation until we find back-scattering form the upper left and right 
corners of our finite simulation box. In this way, we guarantee the correct direction of the 
current along the incompressible strips (see the discussion above). The conductance as a 
function of the magnetic field, the channel distance, and the channel width is estimated 
by calculating the probability N r to find the electron in the upper-right corner within the 
restricted simulation time (see Ref. [2] for details). We made the time propagation by 
using the fourth order Taylor expansion of time evolution operator to the single-electron 
wave function in a 2D real-space grid with the octopus code package |29j . 




Figure 4. Estimated conductance N r as a function of the magnetic flux in the model 
interferometer with different distances between the channels. 



We emphasize that in contrast to the real interferometer (see Fig. [2]), where 
both the channel width and their mutual distance (i.e., the spatial distribution of the 
incompressible strips) is determined by the magnetic field itself, we let in our model the 
magnetic field affect only the flux but not the distribution of the channels. However, for 
each set of calculations, respectively, the width and distance are changed through the 
model parameters given just above. 

In our transport simulations we monitor electron dynamics as a function of 
a magnetic flux added to the background magnetic field (the latter generating the 
incompressible strips). In Fig. |3]we show snapshots of the electron density at different 
times. Here we apply relative fluxes $/$o = and 0.5, respectively. In the beginning 
(see the uppermost row) differences in density between these two cases are not visible, 
so only the other case is plotted. Later on, however, we find a drastic difference between 
the two cases: whereas zero flux leads to transport to the right (second row), a flux 
corresponding to half a flux quantum yields strong current to the left (third row). Thus, 
by exploiting the AB effect, our device is almost completely tunable with respect to the 
current direction at the ouput (left or right). 

To analyze the transport quantization in more detail, we plot in Fig. H] the 
conductance - corresponding here to the electron density transferred to the top of the 
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Figure 5. Estimated conductance N r at zero magnetic field as a function of the 
channel width parameter c for different distances (d) between the channels. 

right channel - as a function of the magnetic flux, and for different distances between 
the channels. We find clear and smooth AB oscillations having exactly the expected 
periodicity. It is noteworthy that in all cases the interchannel transport is completely 
blocked at $/$o — k/2, where k is an odd integer. 

An interesting feature in Fig. H] is the fact that the channel distance has an optimal 
distance d ~ 0.6 when the AB oscillation amplitude is maximized. At smaller (solid 
lines) and larger (dashed lines) values of d the amplitude is decreased. We find that 
the optimal distance corresponds to a case when approximately a half of the density is 
transferred to the right channel at the first (lower) intersection, i.e., the partitioning is 
equal. Therefore, at the second (upper) intersection we find clear AB interference due 
to the phase difference of the optimally partitioned electrons that enclose the given flux. 
Our results agree with the behavior in the AB oscillation strength observed in real AB 
interferometers [30j|3l]- We note, however, that as demonstrated in Fig. [5], the optimal 
distace depends on the channel width c in a nontrivial fashion - apart from the large- 
c limit showing clear suppression of the amplitude. However, it seems that half-half 
partitioning of the density distribution at the first (bottom) point of encounter always 
leads to a relatively high amplitude of the AB oscillation. A complete quantitative 
comparison with experiments would require an accurate determination of the system 
geometry function of the magnetic field. 
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4. Summary 

To summarize, we have performed static and dynamical simulations on Aharanov-Bohm 
interferometers starting from real device parameters. First, the electron density and the 
spatial distribution of the incompressible strips have been obtained self-consistently 
within the Thomas-Fermi approximation. These calculations already suggest that 
interference can take place only if the incompressible strips become leaky, namely not 
strictly compressible or incompressible, and come close to each other, so that partitioning 
of the electron current can take place. These phenomenological considerations have 
then been qualitatively confirmed in the second part of the study, where we time- 
propagate an electronic wave function injected in the channel with tunable parameters. 
We observe distinct Aharonov-Bohm oscillations, whose amplitude strongly depends 
on both the mutual distance between the transport channels and their width. In 
particular, there is an optimal distance yielding maximum oscillation amplitudes. At 
magnetic fields corresponding to half-integer flux quanta the suppression of interchannel 
transport is complete. Taken together, we are able to provide a proof-of-concept for the 
determination of interference patterns in realistic Aharonov-Bohm interferometers in 
the quantum Hall regime. 
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